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Abstract 

The dynamical facilitation scenario, by which localized relaxation events promote nearby relaxation events 
in an avalanching process, has been suggested as the key mechanism connecting the microscopic and the 
macroscopic dynamics of structural glasses. Here we investigate the statistical features of this process via the 
numerical simulation of a model structural glass. First we show that the relaxation dynamics of the system 
occurs through particle jumps that are irreversible, and that cannot be decomposed in smaller irreversible 
events. Then we show that each jump does actually trigger an avalanche. The characteristic of this avalanche 
change on cooling, suggesting that the relaxation dynamics crossovers from a noise dominated regime where 
jumps do not trigger other relaxation events, to a regime dominated by the facilitation process, where a jump 
trigger more relaxation events. 
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1 Introduction 

Structural glasses, which are amorphous solids ob¬ 
tained by cooling liquids below their melting tem¬ 
perature avoiding crystallization, provide an array 
of questions that has been challenging researchers 
in the last decade&ii^i^. These include the nature 
of the glass transition, the origin of the extraordi¬ 
narily sensitivity of the relaxation time on temper¬ 
ature, the Boson-peak, the relaxation dynamics. In 
this respect, here we consider that there is not yet 
an established connection between the short time 
single particle motion, and the overall macroscopic 
dynamics. When observed at the scale of a sin¬ 
gle particle, the motion of structural glasses is well 
known to be intermittent. This is commonly ratio¬ 
nalized considering each particle to rattle in the cage 
formed by its neighbors, until it jumps to a different 
cag©^. Conversely, when the motion is observed at 
the macroscale, a spatio-temporal correlated dynam¬ 
ics emerges^. Dynamical facilitation— by which 
a local relaxation event facilitates the occurrence 
of relaxation events in its proximity, has been sug¬ 
gested as a key mechanism connecting the micro¬ 
scopic and the macroscopic dynamics. Indeed, kinet- 
ically constrained lattice model—, which provide the 
conceptual framework of the dynamical facilitation 
scenario, reproduce much of the glassy phenomenol¬ 
ogy and are at the basis of a purely dynamical inter¬ 
pretation of the glass transition. Different numeri¬ 
cal approaches have tried to identify irreversible re¬ 
laxation events i0iJ^i-J^2ii3a4ii5a6q7p8q9^20 ^ both 

numerical— >2^ and experimental work s revealed 
signatures of a dynamical facilitation scenario. 

Here we provide novel insights into the dynamical 
facilitation mechanisms through the numerical inves¬ 
tigation of a model glass former. We show that it is 
possible to identify single particle jumps that are e/- 
ementar?/relaxations, being short-lasting irreversible 
events that cannot be decomposed in a sequence 
of smaller irreversible events. We then clarify that 
these jumps lead to spatio-temporal correlations as 
each jump triggers subsequent jumps in an avalanch¬ 


ing process. The statistical features of the avalanches 
changes on cooling. Around the temperature where 
the Stokes-Einstein relation first breaks down, the 
dynamics shows a crossover from a high tempera¬ 
ture regime, in which the avalanches do not spread 
and the dynamics is dominated by thermal noise, to 
a low temperature regime, where the avalanches per¬ 
colate. These results suggest to interpret dynamical 
facilitation as a spreading process^, and might open 
the way to the developing of dynamical probabilistic 
models to describe the relaxation of glass formers. 

2 Methods 

We have performed NVT molecular dynamics sim¬ 
ulations^ of a two-dimensional 50:50 binary mix¬ 
ture of 2N = 10^ of disks, with a diameter ratio 
= 1-4, known to inhibit crystallization, at 
a fixed area fraction (f) = 1 in a box of side L. 
Particles interact via an soft potential—, V{rij) = 
e{{aij-rij)/aLT Q{aij - Vij), with a = 2 (Har¬ 
monic). Here rij is the interparticle separation and 
the average diameter of the interacting parti¬ 
cles. This interaction and its variants (character¬ 
ized by different values of a) are largely used to 
model dense colloidal systems, such as foams^, mi- 
crogels^ and glasse s Units are reduced so that 
(Tx, = m = e = fee = 1, where m is the mass of both 
particle species and ks the Boltzmann’s constant. 
The two species behave in a qualitatively analogous 
way, and all data presented here refer to the smallest 
component. 

Cage-jump detection algorithm. We segment the tra¬ 
jectory of each particle in a series of cages interrupted 
by jumps using the algorithm of Ref.—, following 
earlier approaches^. Briefly, we consider that, on a 
timescale <5 of few particle collisions, the fluctuation 
5^(f) of a caged particle position is of the order of 
the Debye-Waller factor (DWF) (rt^). By compar¬ 
ing S‘^{t) with (ti^) we therefore consider a particle 
as caged if S‘^{t) < (m^), and as jumping otherwise. 
Practically, we compute S‘^{t) as {{r{t) — {r{t))s)‘^) 5 , 
where the averages are computed in the time interval 
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[t — 6 : t + 6], with 6 ~ lOt^, and tj, is the ballistic 
time. At each temperature DWF is dehned accord¬ 
ing to Ref.—, (tt^) = {r‘^{tDw))-, where tow is the 
time of minimal diffusivity of the system, i.e. the 
time at which the derivative of log(r^(f)) with re¬ 
spect to log(t) is minimal. At each instant the al¬ 
gorithm allows to identify the jumping particles and 
the caged ones. We stress that in this approach a 
jump is a process with a finite duration. Indeed, by 
monitoring when 5^ equals (u^), we are able to iden¬ 
tify the time at which each jump (or cage) starts and 
ends. We thus have access to the time, tp, a parti¬ 
cle persists in its cage before making the first after 
an arbitrary chosen t = 0 (persistence time), to the 
waiting time between subsequent jump of the same 
particle (cage duration), and to the duration Atj 
and the length Arj of each jump. 

3 Results 

3.1 Jumps as irreversible elementary pro¬ 
cesses 

The idea of describing the relaxation of structural 
glasses as consisting of a sequence of irreversible 
processes is not new, and different approaches have 
been followed to identify these events. For instance, 
irreversible events have been associated to change 
of neighbor s ^ displacements overcoming a 

threshold in a fixed time laps^, to processes iden¬ 
tified through clustering algorithm applied to the 
particle trajectorie s or to more sophisticated 
approache&l^. We notice that since at long time 
particles move diffusively, all procedures that coarse 
grains the particle trajectory enough will eventually 
identify irreversible events. Here we show that the 
jumps we have identified are irreversible, and we give 
evidence suggesting that these can be considered as 
‘elementary’ irreversible events, i.e that they are the 
smallest irreversible single-particle move, at least in 
the range of parameters we have investigated. 

Investigating both the model considered her©^, 
as well as the 3d Kob-Andersen Lennard-Jones (3d 



Figure 1: Average persistence time, (tp), cage du¬ 
ration, (tw) and jump duration, (Atj), as a func¬ 
tion of the temperature, (t^) grows as an Arrhe¬ 
nius (tw) oc exp (A/T) (red full line), whereas {tp) is 
compatible with several super-Arrhenius laws. The 
black full line is, for example, a fit {tp) oc exp {A/T'^), 
while the black dashed line is a Vogel-Fulcher law 
{tp) oc exp {B/{T — Tq)), predicting a divergence at 
a finite temperature Tq ~ 0.001. The arrow indicates 
the temperature Tj, = 0.002 where {tp) and {t^)) de¬ 
couple and the SE relation breaks down. Conversely, 
{Atj) remains roughly constant on cooling. 

KA LJ) binary mixture^ and experimental colloidal 
glass^, we have previously shown that the protocol 
dehned in Sec. [2] leads to the identihcation of irre¬ 
versible events. Indeed, the mean square displace¬ 
ment of the particles increases linearly with the num¬ 
ber of jumps, allowing to describe the dynamics as a 
continuous time random walk (CTRW)— . 

Within this approach two fundamental timescales 
are found, the average persistence time {tp) and the 
average cage duration {t^)- The former corresponds 
to the relaxation time at the wavelength of the order 
of the jump length (Arj), while the latter is related 
to the self diffusion constant, D oc {Ar'j)/{t^). Fig{T] 
shows that the two timescales are equal at high tem¬ 
perature, but decouple at a temperature 0.002, 
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Figure 2: Mean squared jump length {A.rj) as a func¬ 
tion of the jump duration Atj at different tempera¬ 
tures. 

which marks the onset of the Stokes-Einstein (SE) 
breakdown at the wavelength of the jump length. We 
find that shows an Arrhenius temperature de¬ 
pendence (t^) (X exp (A/T), while (tp) increases with 
a faster super-Arrhenius behaviour (see the caption 
of Fig{T]). It is worth noticing that the decoupling 
between the average persistence and waiting time, is 
known to control the breakdown of the SE relation 
at generic wavelengths, and to induce temporal het¬ 
erogeneities^*^. These Endings suggest that Tx may 
represent a crossover from a localized to a more cor¬ 
related relaxation process. A similar scenario has 
been recently reported for models of atomic glass 
forming liquids, where the SE breaks down and the 
size of dynamics heterogeneities markedly accelerates 
below a well defined value of Tx-— 

We performed two investigations supporting the 
elementary nature of the jumps we have identified. 
Eirst, we have considered the change of the average 
jump duration (Atj) on cooling, as the duration of 
elementary relaxations is expected not to grow with 
the relaxation time. Eig. [T] shows that the (Atj) 
is essentially constant, despite the relaxation time 
(tp) varying by order of magnitudes. Indeed, at low 
temperature {tp)/{Atj) ^ 1, clarifying why we call 
them ‘jumps’. Then we have considered how parti- 



Eigure 3: Excess probability to observe contempo¬ 
rary jumps, Cj{r, 0), as function of the distance and 
at different temperature, as indicated. The dashed 
line is a guide to the eyes oc exp(—1.35/r). 

cles move while making a jump. Fig. [2] illustrates 
that the mean squared jump length grows subdif- 
fusively as a function of the jump duration, with 
a subdiffusive exponent that decreases on cooling. 
Conversely, one would expect a diffusive behaviour 
if jumps were decomposable in a series of irreversible 
steps. 

These results supports the identification of the 
jumps we have defined with the elementary relax¬ 
ations leading to the macroscopic relaxation of the 
particle system. 

3.2 Correlations between jumps 

While each particle behaves as a random walker as 
it performs subsequent jumps, yet jumps of differ¬ 
ent particles could be spatially and temporally corre¬ 
lated. We investigate these correlations focusing on 
the properties of a jump birth scalar field, dehned as 

1 ^ 

=—'^bi{t)5{r-nit)). (1) 

i 

Here bi{t) = 1 if particle i starts a jump between t 
and t+5t, where 5t is our temporal resolution, bi[t) = 
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0 otherwise. The scalar field b allows to investigate 
the statistical features of the facilitation process by 
which a jump triggers subsequent ones. To this end, 
we indicate with (6(r, t))6(o,o)=i the probability that 
a jump starts in {t, r) given a jump in (t = 0, r = 0), 
and investigate the correlation function 


Cj{r,t) 


- ( 6 ) 

9{r, t) 


( 2 ) 


Here g{r, t) is a time dependent generalization of the 
radial distribution function 

= 27rrp(N - 1) ^ ~ ~ 

(3) 

through which we avoid the appearance of spurious 
oscillations in the correlation function Cj (r, t) due 
to the short range ordering of the system. In EqjSj 
(6) is the spatio-temporal average of the jump birth, 
and decreases on cooling as (6) = + (Atj))“^ 

(at low temperature {b) ~ {tiu)~^ as << {Atj)). 
Accordingly, the correlation function Cj{r,t) is the 
probability that a jump triggers a subsequent one at 
a distance r after a time t. 

We first consider the spatial correlations between 
contemporary jumps, where two jumps are consid¬ 
ered contemporary if occurring within our tempo¬ 
ral resolution. Fig. [3] shows that C'j(r, 0) decays ex¬ 
ponentially, with a temperature independent corre¬ 
lation length ^j(0,T) ~ 1.35. This result clarifies 
that jumps aggregate in cluster of roughly Worr — 
p7r^j{0) ~ 5 events. A similar scenario has been ob¬ 
served in a different model system, where jumps have 
been observed to aggregate in clusters of roughly 7.6 
particles^. Our results also support previous find¬ 
ings suggesting— that the elementary excitations of 
structural glasses have a temperature-independent 
length not larger than few particle diameters and are 
consistent with a recently introduced first principle 
extension of the Mode Coupling Theory—. The in¬ 
vestigation of the displacements of the particle jump¬ 
ing in each cluster does not reveal characteristic spa¬ 
tial features. Structured particle motion, such as 
string-like particle displacements—^ or displacements 



Figure 4: Evolution of the spatial correlation be¬ 
tween jumps with time. Each panel refer to a 
different temperature, as indicated. Within each 
panel, the different curves correspond to t = 
0,10,20,30,100,500,10^ and 10^, from top to bot¬ 
tom. At high temperature data corresponding to 
the largest times are missing as the correlation is 
too small to be measured. 


reminiscent of T1 events^ must therefore result from 
a succession of events rather than a single one. 

We now consider the time evolution of the spatial 
correlation between jumps. Fig. [5] illustrates that 
at all temperatures and times the decay of the cor¬ 
relation function is compatible with an exponential, 
Cj{r,t) (X A(t) ex.p{—r/^j{t)). The time dependence 
of the amplitude is illustrated in Fig. [5l At all tem¬ 
peratures the short time decay of the amplitude is 
exponential, A{t,T) = A(0, T) exp(—t/ryi(T)), the 
characteristic decay time slightly increasing on cool¬ 
ing. While no other decay is observed at high tem¬ 
peratures, at low temperatures the exponential decay 
crossovers towards a much slower power-law decay 
A{t) ~ with a ~ 0.4. Fig. [6] shows that the cor¬ 
relation length slowly grows in time, approximately 
as ~ with b ~ 0.1. 

The initial fast decrease of the amplitude makes 
difficult to obtain reliable estimates of its time de- 
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Figure 5: Panel a shows that the amplitude A{t) of 
the jump correlation function Cj{r,t). Panels b and 
c clarify that a hrst exponential decay is followed, 
at low temperature, by a second one, which approx¬ 
imately follows a power law. 



Figure 6: Time dependence of the jump correlation 
length, at different temperatures. The data suggest 
that at low temperature the correlation length slowly 
grows in time, as ij{T) oc 


pendence and correlation length, despite intense 
computational efforts. Nevertheless, our data 
clearly show the reported exponential to power-law 
crossover in the decay of the amplitude of Cj{r,t). 


The highest temperature at which this decay ex¬ 
hibits a power law tail, is consistent with the tem¬ 
perature Tx where {tw) and {tp) first decouple, and 
the SE relation breaks down (see Sec J3.1l) . This 
suggests that the breakdown of the SE relation is 
related to a crossover in the features of the facil¬ 
itation process. We investigate this crossover fo¬ 
cussing on the number of jumps triggered by a given 
jump. This is given by Ntr{T) oc n(t,T)dt, 
where n{t,T) = f C(r,t)rdr oc A{t,T)^‘^(t,T)dt, 
is the number of jumps it triggers at time t. As 
at high temperature the variation of the correlation 
length is small with respect to that of the ampli¬ 
tude, one can assume — C(0,T) and esti¬ 

mate Ntr{T) oc A(0,r)^^(0,T)r^(r). At low tem¬ 
perature, the integral is dominated by the long time 
power law behavior of the amplitude and of the cor¬ 
relation length, and the number of triggered events 
diverges as Atr(E, t) oc A{t)^‘^{t)dt oc ^-“+2^+1 qc 
fO-S. 

4 Discussion 

We conclude by noticing that the above scenario sug¬ 
gests to interpret facilitation as an infection spread¬ 
ing process, in which a particle is infected each 
time it jumps. Since each particle can be infected 
more than once, the relevant infection model is of 
susceptible-infected-susceptible (SIS) type. In this 
framework, the exponential to power-law crossover 
in the decay of the amplitude of Cj{r,t) signals a 
transition from a high temperature resilient regime, 
in which a single infected site only triggers a finite 
number of infections, to a low temperature regime 
in which the number of triggered infection diverges. 
A complementary interpretation can be inspired by 
the diffusing defect paradigm— ii. We suggest that 
the correlation length of contemporary jumps, 0(0); 
is akin to the typical defect size, which, according to 
our results, is temperature independent. In the high 
temperature regime, this is the only relevant corre¬ 
lation length, as defects are rapidly created and de¬ 
stroyed by noisy random fluctuations, before they 


6 


















can sensibly diffuse. At low temperature, the ef¬ 
fect of noise becomes smaller: the short time cor¬ 
relation length is still dominated by the defect size, 
< ta) — 0(0)) whereas its long time behaviour, 
0(i >> ta), is controlled by the typical distance 
defects have moved up to time t. Further studies are 
necessary to investigate which of the two interpreta¬ 
tions is more appropriate. 
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